Measurement and simulation of laser-induced fluorescence from non-equilibrium 

ultracold neutral plasmas 
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We report new measurements and simulations of laser-induced fluorescence in ultracold neutral 
plasmas. We focus on the earliest times, when the plasma equilibrium is evolving and before the 
plasma expands. In the simulation, the ions interact via the Yukawa potential in a small cell with 
wrapped boundary conditions. We solve the optical Bloch equation for each ion in the simulation 
as a function of time. Both the simulation and experiment show the initial Bloch vector rotation, 
disorder-induced heating, and coherent oscillation of the rms ion velocity. Detailed modeling of 
the fluorescence signal makes it possible to use fluorescence spectroscopy to probe ion dynamics in 
ultracold and strongly coupled plasmas. 

PACS numbers: 32.50.-|-d, 52.20.-j, 52.25.0s, 52.65.Rr, 52.27.Cm 



I. INTRODUCTION 



Ultracold neutral plasmas can be formed by photo- 
ionizing laser-cooled gases P, 0, H, S S @| • The initial 
ion temperature is approximately equal to the the laser- 
cooled atom temperature, typically 0.01 K or less. These 
plasmas are also small. Their characteristic size is typ- 
ically 1 mm. The electron temperature can be tuned 
from approximately 0.5 K to 1000 K. At low initial elec- 
tron temperature, space-charge effects prevent electrons 
from leaving, and the plasmas are charge neutral [l| . The 
electrons screen the ion-ion interactions, playing an im- 
portant role in the rate at which the ions thermalize the 
the plasma _eyolvesJ^J^ @J^Jl^Jll|Jl|Jl^^ M, M, 
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Immediately after photo-ionization, the ion-ion inter- 
action raises the ion temperature. This happens because 
the ions are randomly distributed throughout the plasma. 
Although they have essentially zero kinetic energy, they 
have a great deal of electrical potential energy. The ions 
move to minimize this potential energy, a process called 
disorder- induced heating [1, . The initial ion motion 
is coherent in a sense, and damped oscillations have been 
observed in the ion temperature evolution. 

During this initial plasma phase, the ion temperature 
is not well defined [20|. This is somewhat surprising 
because the initial nearest-neighbor distribution, which 
dominates the disorder induced heating phase, is approxi- 
mately Gaussian. It has also been shown that the spatial 
distribution evolves self-similarly (i.e.: is always Gaus- 
sian if the ion temperature is not too low) [T^j- Further- 
more, the velocity distributions both before and after the 
disorder-induced heating phase are Gaussian. Therefore 
one might expect a nominally Gaussian distribution of 
velocities throughout the thermalization time. 

Early fluorescence and absorption studies of ions in 
these plasmas used a Voigt profile to extract the ion tem- 



perature 0,13. The Voigt profile, 

S'(i^) = / L[v — x )G{x') dx , 



(1) 



is a convolution of the Lorcntzian natural line shape L{i') 
with the Gaussian G{v) thermal line shape. While this 
appears to be valid in the later stages of the plasma evo- 
lution, the Gaussian width has limited meaning as the 
ions thermalize. Recent work has shown that Eq. [T]may 
be valid locally in the plasma [28j . 

In this paper, we present a study of laser-induced flu- 
orescence in ultracold neutral plasmas. Rather than ex- 
tract the temperature using a lineshape analysis or local 
fluorescence imaging, we perform a simulation of the flu- 
orescence by numerically integrating the optical Bloch 
equations for the plasma ions. The simulated fluores- 
cence signal is compared with experimental results. Both 
the experiment and simulation show the initial rotation 
of the Bloch vector as the ions (all initially in the ground 
state) begin to absorb and scatter laser light. During 
the disorder-induced heating phase, the linewidth broad- 
ens and the fluorescence signal falls. At high densities, a 
coherent oscillation in the ion motion is observed. The 
experiment and simulation have satisfactory numerical 
agreement, and we suggest ways to use these compu- 
tational techniques in future experiments in strongly- 
coupled neutral plasmas. 



II. EXPERIMENT 

Up to seven million "^^Ca atoms are laser-cooled and 
trapped in a magneto-optical trap (MOT). The trap den- 
sity is approximately Gaussian of the form 



n{r) — no exp (— r^/2(T^) , 



(2) 
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where uq is the peak density with a typical value of 10 
cm~'^ in our experiment, and a is the one-dimensional 
rms size of the MOT, typically 350 ^m. 
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FIG. 1: Partial energy level diagram, transition wavelengths, 
and transition probabilities for relevant levels in (a) Ca I and 
(b) Ca II. The initial electron energy Ee is determined by 
the energy difference between the photoionization laser pho- 
ton energy and the Ca I ionization energy (I. P.). The initial 
electron energy is Te = 2i?e/3fcs. The 423 nm Ca I transition 
is also used for the calcium MOT. Not shown is a repumper 
laser at 672 nm that plugs an optical leak from the MOT. 
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FIG. 2: Relative fluorescence signals from a high density 



(black line, no = 7.3 x 10 cm~ ) and low density (gray line, 
no = 0.45 X 10^ cm~^) plasma at Te = 57 K. In the high den- 
sity plasma, disorder induced heating occurs more rapidly. 
The fluorescence signal shows a weak oscillation around 100 
ns. At 500 ns the signal has decayed somewhat less completely 
compared to the low density measurements. All fluorescence 
curves are in plotted in arbitrary units (linear scale) and have 
been normalized with a peak of 1. 



Atoms in the trap are photoionized using ns-duration 
pulsed lasers at 423 and 390 nm. The temperature of 
the plasma ions is approximately equal to the ^mK tem- 
perature of the MOT atoms. The initial energy of the 
plasma electrons is determined by the wavelength of the 
390 nm laser. It is equal to the difference between the 
laser photon energy and the atomic ionization potential 
(see Fig. [T]). 

After the plasma is generated, plasma ions are excited 
using a cw laser beam tuned to the 4s — > 4p 
transition at 397 nm. Fluorescence at this wavelength is 
collected using a lens, isolated using a band-pass optical 
interference filter, detected using a photomultiplier tube, 
and recorded using a fast digital oscilloscope. While most 
of the fluorescence from the upper state is at 397 

nm, there is also a 6% branch to a dark, metastable 
3d state. The 397 nm laser beam is focused to a 

Gaussian waist of 100 /im in the center of the plasma. 
Typical fluorescence curves are shown in Fig. [2l 

Because the ions are all initially in the ground state, 
it takes time for them to start scattering photons from 
the probe laser beam. In the fluorescence signal, this is 
seen as the rapid initial rise at time i = 0. As mentioned 
already, the ions initially accelerate due to the forces from 
other ions in the plasma. This leads to disorder-induced 
heating. This causes a rapid fall in the fluorescence signal 
after about 30 ns. 

The plasma also expands because it is no longer con- 
flned by the MOT. When the initial electron energy is not 
too low, this plasma expansion velocity is approximated 



fairly well by the equation 

Wcxpl?', t)^r — t, (3) 

where r is the radial coordinate, /cg is Boltzmann's con- 
stant, Te is the electron temperature, m is the ion mass, 
and t is time. This equation is valid only for very 
early times in the plasma evolution, before the plasma 
has expanded, and before the electron temperature has 
changed due to expansion, evaporation, or recombina- 
tion. This expression also neglects the influence of cor- 
relations, velocity-changing collisions, or many-body ef- 
fects. A characteristic time for plasma expansion is when 
I'cxp^ = c. Substituting the rms value of r = v^c, we 
find that for Te = 57 K, the characteristic time is approx- 
imately 2 /is. As the plasma expands, the fluorescence 
continues to fall. In our simulation, we can neglect the 
overall plasma expansion by concentrating on the flrst 
few hundred ns of data. 

As can be seen in Fig. [2l the fluorescence method 
clearly distinguishes between high and low density. 
Disorder-induced heating occurs more rapidly in high 
density plasmas than in low density plasmas. The flu- 
orescence signal shows a weak oscillation around 100 ns 
in high density plasmas that is absent in low density plas- 
mas. Finally, at 500 ns, the fluorescence signal does not 
decay as completely for high density plasmas as it does for 
low density plasmas. For a fixed density, the fluorescence 
signal also changes as the initial electro temperature is 
reduced. If this signal can be properly interpreted, de- 
tails about electron shielding, three-body recombination, 
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many-body effects, and particle correlations can be de- 
termined at very early times in the plasma evolution. 



The fraction of atoms in the excited state is given by ■ 
For example, if there is no 397 nm laser, then = and 
the equation for is 



III. SIMULATION 



(i) = [1 



(0)] 



1. 



(11) 



We simulate the plasma ion fluorescence signal by in- 
tegrating the optical Bloch equations for a collection of 
ions in a cell. In the simulation, ions interact via the 
Yukawa potential, 



V{t) = 



Alien 



(4) 



where go is the fundamental charge, Aj^ = ^JkBT^eQ/nq^ 
is the Debye length, fcs is Boltzmann's constant, Te is 
the electron temperature, n is the density, and eq is the 
permittivity of free space. 

Plasma ions are randomly distributed over a cubic cell. 
The cell dimension is L « cr/ 10. For the Gaussian den- 
sity profile of Eq. [21 the density is constant across this 
cell size. We differentiate Eq. |3]to find the a;,?/, and z 
components of the force on each ion due to the presence 
of all of the other (screened) ions in the cell. We use a 
fourth-order Runge Kutta stepper to move the ions in 
time. 

At each time step, we also solve the optical Bloch equa- 
tions for each ion in the cell. The ions are approximated 
as two-level atoms. The Hamiltonian for an ion is written 
in terms of the detuning lo and the Rabi frequency f2 as 



— [a- + a+) 



where the the raising and lowering matrices are 
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(5) 



(6) 



The equation for the density matrix is 

/9 = — iLO [a +(J ^ p — pa -\-a ^) 

- i^[{a++a-)p- p{a++a^)] (7) 

+ 2 (^^-^'^+ ^'^+^-^^^'^+'^-)■ 
where 7 is the total decay rate of the 4p state. 
This equation of motion is somewhat more compactly 
written using the Pauli spin matrices, with cfx ~ CT+-l-cr_, 
Gy = — i(T+ + i(T_, and ct^ the 2x2 identity matrix: 



dt 
d{<Jv) 

dt 
d{cjx) 

dt 



= (a,) -7(1 + (a,)) 

= UJ {(Jx) - n {(Jz) - ^ {(Ty) 



(8) 
(9) 
(10) 



With the 397 nm laser present, and the fluores- 
cence signal depends on time as 



/W = i[l 



(12) 



Because the Bloch equations depend on the detuning uj 
of the laser beam from the 397 nm resonance transition, 
we need to consider two frequency shifts. The first is 
the motion of the ions in the cell. This is given by the 
first-order Doppler shift. 



2ttu 



(13) 



where u is atomic transition frequency in the rest frame of 
the atom, c is the speed of light, and v is the component 
of the atomic velocity along the direction of the laser 
beam propagation. 

The second frequency shift comes from the radial ac- 
celeration of the plasma (Eq. [3]) . For early enough times 
in the plasma evolution, the ions accelerate and begin to 
move, but the overall density profile is unchanged. There- 
fore we can simulate the frequency shift due to plasma 
expansion by placing our cell at different plasma radii, 
calculating the mean velocity according to Eq. [31 and 
adding this directed velocity to all of the ion velocities in 
the cell. 

For each time step, we calculate the position and ve- 
locity of each ion. We use wrapped boundary conditions 
to maintain a constant number of ions in the cell. We 
use the velocity (including the effect of plasma expan- 
sion) to determine the ion's frequency shift uj and use 
a second-order Runge-Kutta routine to calculate its (Tz- 
The (Tz is then averaged over the ions in the cell to simu- 
late the fluorescence signal. The simulation is performed 
for many different initial ion densities. The ion density 
is chosen randomly from a Gaussian distribution so that 
the different density regions are appropriately weighted. 

Care must be taken to properly simulate the fluores- 
cence decay into the optically dark 3d state. At each 
time step the total fluorescence rate from the level 
is calculated and multiplied by the branching ratio to get 
the decay rate into the dark 3d ^-D3/2 level. This decay 
rate is multiplied by the time step dt to get the prob- 
ability that the ion has made a transition to this level. 
As is traditionally done in Monte Carlo simulations, this 
probability is compared to a random number with a flat 
distribution between and 1. If the random number is 
less than the probability, then the simulated ion makes 
the transition to the dark state and no longer fluoresces. 
Decay into the dark metastable state reduces the fluores- 
cence signal, particularly at longer times. 
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FIG. 3: Comparison of experiment (solid line) with simulation 
(dashed line) for a range of initial plasma densities. Both the 
simulation and experiment show initial rotation of the Bloch 
vector into the excited state, disorder-induced heating, and 
plasma expansion. As the density increases, the simulation 
shows the oscillation in the rms velocity distribution occur- 
ring later in time than in the experiment. This may be due 
to the onset of many-body interactions. For all of these mea- 
surements, the initial electron temperature is Te = 57 K. 



IV. COMPARISON WITH EXPERIMENT 

The simulation is compared with the experimental flu- 
orescence signal in Fig. |3l The three phases of the signal 
are clearly visible. The first is the rapid rise of the fluo- 
rescence. This corresponds to the initial rotation of the 
z-component of the Bloch vector into the excited state. 
The second is an almost equally rapid decline in fluores- 
cence after 20 ns. This is due to disorder-induced heating 
- the acceleration of the ions as they move to minimize 
their electrical potential energy. The third is the slower 
signal decay at longer times. This is due to the overall 
plasma expansion and to shelving into the Zd state. 

The agreement between the simulation and experiment 
is quite good. It suggests that the essential physics of the 



ion acceleration and plasma expansion are captured in 
this relatively simple approach. In particular, the initial 
Bloch vector rotation and disorder-induced heating are 
well-simulated. 

There is an apparent discrepancy in the ion tempera- 
ture oscillation at high density. As has been noted pre- 
viously, the ions are all initially at rest. They begin to 
move in their local electrical potential energy well. After 
a time approximately equal to the inverse of the plasma 
frequency, the ions on average have moved to the other 
side of their local well and come to rest as they "collide" 
with their nearest neighbor. This motion - initially at 
rest, simultaneous acceleration, coming to rest again af- 
ter one plasma period - produces an oscillation in the rms 
ion velocity distribution. Because nearest neighbors are 
randomly located in the initial plasma, the local poten- 
tial is not parabolic and the collision time is not exactly 
the same for the ions. After a short time, the coherence 
in the rms ion velocity oscillation is lost and the fluores- 
cence curve smooths out. 

This discrepancy may be due to collisions or opti- 
cal pumping in the plasma. Velocity changing colli- 
sions and additional shifts from the MOT magnetic fields, 
and many-body interactions may also influence the flu- 
orescence signal. These possibilities are currently under 
study. 



V. CONCLUSION 

We present measurements and simulations of laser- 
induced fluorescence of ions in an ultracold neutral 
plasma. In the simulation, the ions interact via the 
Yukawa potential. The ion position and velocity is calcu- 
lated as a function of time, including the influence of ions 
in the cell as well as the overall plasma expansion. At 
each time step we solve the optical Bloch equations and 
track the excited-state population to simulate the fluores- 
cence signal. We find good agreement between simulation 
and experiment. 

Fluorescence measurements can be used to probe the 
early-time ion dynamics over a wide range of initial 
plasma densities. The ability to accurately calculate fluo- 
rescence properties of ultracold plasmas makes it possible 
to quantitatively interpret laser-induced fluorescence sig- 
nals from ultracold plasmas. Exploring the influence of 
an external magnetic field, probe laser detuning, varying 
optical depth in the plasma, displacing the probe laser 
beam from the center of the plasma, and other impor- 
tant effects can be done as straightforward extensions of 
our approach. 

Future work will focus on the agreement at higher den- 
sities in the plasma. This regime is important to under- 
stand because it is where many-body effects are likely 
to appear. Similarly, work at lower initial electron tem- 
perature is also interesting for the same reason. At the 
very lowest initial electron temperature or at the high- 
est plasma densities, the energy scale in the simulation 
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is determined only by the disorder-induced heating, both 
for electrons and ions. Comparing the Yukawa simula- 
tion, the experiment, and a direct molecular dynamics 
calculation would provide good insight into the onset of 
correlations in this regime. 
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